Electroproduction of the d* dibaryon 
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The unpolarized cross section for the electroproduction of the isoscalar = 3^ di-delta dibaryon d* 
is calculated for deuteron target using a simple picture of elastic electron-baryon scattering from the 
AA(^Z)i) and the NN{^Si) components of the deuteron. The calculated differential cross section 
at the electron lab energy of 1 GeV has the value of about 0.24 (0.05) nb/sr at the lab angle of 
10° (30°) for the Bonn B potential when the dibaryon mass is taken to be 2.1 GeV. The cross 
section decreases rapidly with increasing dibaryon mass. A large calculated width of 40 MeV for 
d*(AA ^Sa) combined with a small experimental upper bound of 0.08 MeV for the d* decay width 
appears to have excluded any low-mass d* model containing a significant admixture of the AA(^S'3) 
configuration. 



PACS numbers: 14.20.Pt, 25.30. Dh, 13.85.Fb 
I. INTRODUCTION 

Given the availability of excellent electron beams, the 
electroproduction of dibaryon resonances represents a 
promising way to look for these long-soug ht objects §,|. 
In assessing the practical prospects of such experiments, 
one needs to begin with rough estimates of the electro- 
production cross sections. This has to be done separately 
for each dibaryon candidate of interest, since each candi- 
date may involve unique theoretical issues. 

One candidate dibaryon that has been discussed re- 
cently is the di-delta d* of quantum numbers J'^T = 3+0 
iQ-Q. This dibaryon is a six-quark state which at large 
inter-baryon separations may be visualized as a pair of 
A's. Preliminary studies of its electroproduction from 
deuteron targets have recently been made by Qing |^ and 
by Sun Q, both of Nanjing University. They find that 
the KroU-Ruderman AA production process does 

not contribute because of the special isospin structure of 
d and d* . They concentrate instead on another process 
important in the electroproduction of two pions from pro- 
ton |l^ and deuteron |l^ based on the electroexcitation 
of a nucleon isobar N*{1520) in the intermediate state, 
namely jd NN * (1520) AA. This mechanism can 
be interpreted as the production of d* through its NN* 
component. 

However, these preliminary results for the calculated 
differential cross sections turn out to be very small at 
small momentum transfers (say < 500 MeV/c), almost 
two orders of magnitude smaller than a simple estimate 
made by me using a simple picture of inelastic scatter- 
ing from the AA(^Z?i) of the deuteron, as described in 
diagram a of Fig. |l|. These results seem to suggest 
that d* is more easily produced from the deuteron at 
these small momentum transfers through elastic electron- 
baryon (eB) scatterings rather than inelastic excitations 
of the struck baryon. 

It is worth noting in connection with diagram |l^ that 
if the d* is visualized as a ^6*3 didelta, there is no produc- 



tion from the deuteron AA('^S'i) component in the low- 
est order. This is because the electromagnetic operators 
are of rank only or 1 in intrinsic spin. Hence electro- 
production in this simple lowest-order picture is possible 
only from the AA(^Z3i) component of the deuteron. 

In addition to diagram or equivalently its per- 
turbative part diagram ^d, d* can also be electropro- 
duced through the NNi^D^ — G3) component of d*, 
as described in diagram np. Lomon has looked for a d* 
resonance in these channels by studying the energy de- 
pendence of their NN phase parameters . These 



NN{'^D3 — Ga) channels are in fact the dominant de- 



cay channels of d* [^||ljJT^. He has found that these 
phase parameters are consistent with the absence of any 
dibaryon resonance with a width exceeding 1 (2) MeV 
near 2.1 (2.25) GeV. 

A much stricter upper limit on the width of an np 
dibaryon resonance has been obtained experimentally 
some time ago by Lisowski et al. p6[ | . They measured the 
total np cross section at cm. energies 2.00 to 2.23 GcV 
at high energy resolution, namely about 1.4 MeV at 2.11 
GeV. No evidence was seen for narrow resonances with 
areas greater than 5 mb-MeV. At 2.1 GeV where the np 
I — total cross section is 33 mb, a totally elastic "^Z^a 
resonance would have a maximum unitarity-limited cross 
section of crmax = 39 mb. Its energy-integrated area is 
(7r/2)(Tniaxr for a pure Breit-Wigner shape, where T is 
its total decay width. The assumption of a pure Breit- 
Wigner resonance should be a very good one, because 
the NN phase shift at the total energy of 2.1 GeV 
is only 4° ||l|. The experimental bound of 5 mb-MeV 
on the integrated area translates to an upper bound of 
only 0.08 MeV for the width of any resonance that has 
escaped detection in this experiment. 

The theoretical decay width of d* has been estimated 
in 0|. If the d* is taken to be a bound system of two 
finite-sized A's, describable by a nonrelativistic quark 
model, the decay has been estimated to be about 10 
MeV when its mass is 2.1 GeV. The result is, however, 
sensitive to both dynamical input and rescattering cor- 
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rections; it could be as large as 40 MeV for a realistic po- 
tential treated perturbatively without any rescatterinfg 
correction, as I shall show in more detail in this paper. 

It has been argued in Q that if the dibaryon is made up 
of six delocalized quarks instead of two separated 

A's, its decay width could well be smaller by an order 
of magnitude, say 1 MeV. This suggestions remains to 
be confirmed by a detailed calculation. A reduction to 
0.1 MeV would be much harder to realize. Furthermore, 
if the estimated decay width has to be reduced for any 
reason, its estimated production cross sections, including 
the electroproduction cross section estimated here using 
a didelta model, must also be reduced correspondingly. 

Thus the case for a d* resonance near 2.1 GeV does not 
appear promising. However, the d* could have a mass 
higher than expected if the dynamnics is different from 
that described by the delocalization and color-screening 
model. Hence it is still of interest to study the electro- 
production cross section for c?*, should such a resonance 
exist. 

It is obvious that in addition to diagram |l|a or b, one 
should also included the process shown in diagram |l|c, 
where the d* is produced through the NN{'^D3) channel, 
the NN{^G3) contribution having been ignored in this 
lowet-order picture. 

Although these two processes alone do not add up to a 
quantitative description of the electroproduction, they 
are of sufficient interest to justify the detailed report 
given here, including the contributions of convective and 
magnetization currents. The picture is necessarily very 
rough, because of the neglect of many components and 
processes. However, uncertainties about the mass and 
structure of d* and about short-distance nuclear dynam- 
ics have discouraged me from undertaking a more ambi- 
tious calculation at this time. 

The paper is organized as follows: The notation used 
is defined in Sec. II where brief comments relevant to 
the present calculation are made. Sec. Ill shows how 
the calculation is done for diagram |l]a, when the needed 
AA(^_Di) component of the deuteron is already available. 
The contributions from different reduced matrix elements 
(RME's) are briefly discussed. 

Sec. IV shows that the inclusion of the d*{NN^D3) 
contribution via the perturbative process Fig. Wp can be 
broken down into three steps: (a) a perturbative eval- 
uation of the d*{NN^D3) wave function, (b) a calcula- 
tion of the d* decay width for the same input dynamics 
(taken to be the 28-channel Argonne potential jl^] and 
the Bonn B potential |l^ ) , and finally (c) the evaluation 
of the electroproduction cross section itself. 

Calculational results are presented in Sec. V, where 
the production amplitude from diagram ^c is found to 
be greater than that from diagram |^a by a factor of 2- 
4 at certain angles or momentum transfers. At a lab 
angle of 30° and an electron lab energy of 1 GeV, I find a 
differential production cross section of about 0.05 nb/sr 



when integrated over the energy loss. 

My calculated cross section is smaller by two orders 
of magnitude from the preliminary value of about 10 
nb/sr obtained recently by the Nanjing group (to be 
called QSW) ||l^ using a different method of calcula- 
tion. Part of the discrepancy, accounting for a factor 
of 5, comes from the fact that QSW has included only 
TT exchange but not p exchange. The remaining discrep- 
ancy, of more than an order of magnitude, must be due to 
differences in calculational methods used. For example, 
I use a principal-value Green function for two nucleons 
in the d*{NN^D3) wave function, whereas QSW use an 
outgoing-wave boundary condition. To help in disentan- 
gling the discrepancy in the future, I have included many 
details in the present paper. 

The decay width calculated here for the AA(^53) 
model of d* is about 40 MeV at m* =2.1 GeV. When 
used with the experimental upper bound of only 0.08 
MeV of any / = 0, J = 3 up resonance ||l^, the calcu- 
lated width makes it very unlikely that any undetected 
d* has a significant probablity of the AA(^5'3) configu- 
ration. 



II. ELECTROPRODUCTION CROSS SECTION 

FOR ED ^ ED* 

Consider the scattering of a relativistic electron beam 
of lab energy e to the lab angle 6* by a target that goes 
from an initial state z to a final state / as the result of 
the scattering. The differential cross section in the lab 
(after integrating over the energy transfer) is known to 
have the form |2(J 



an 



aMf7,l[vLTt+VT{Tt + Tr^) 



where 



Cm = 



acos{e/2) 



2esin''(6l/2) 



(1) 



(2) 



is the Mott differential cross section, a is the fine struc- 
ture constant. 



/r, 



1 



2esin2(6l/2) 



target 



is a target recoil factor. 



q- 



tan 



(3) 



(4) 



are the electron kinematical factors that depend on its 
four-momentum transfer Q = K—K' = [lo, q). Here K = 
(e, k) is the four- momentum of the electron in the lab 
before the scattering, and K' its four-momentum after 
the scattering. The notation is that of pl[, with < 0. 
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The target factors are 
An 

J=0 



Ti 



2.L 



-J2\{Jf\\Mj{q)U) 



fi — 



47r 



2J,: 



(5) 



J=l 



where a — el or mag. The nuclear reduced matrix ele- 
ments (RME) {Jf\\0\\Ji) that appear are all dimension- 
less quantities. The operators O are the longitudinal 
Coulomb, and the transverse electric and magnetic, mul- 
tipole operators that have the following simple forms in 
momentum space [|| 



Aw 



111 



rr'iq) 



An 
An 



dr!q[i;(17q)® J(q)]'^, 

dn^[Y,in^)^3{q)y . (6) 



Here Yj is a spherical harmonic, p(q) and J(q) are the 
baryon charge and current density operators. 

The current density operator can be separated into 
convective and magnetization terms: 



J(q)-Jc(q)+Js(q). 



(7) 



The spin magnetization term Js = V x /ig originates 
from the dibaryon isoscalar magnetic moment operator 
As = MoMsS. Here /io is the nuclear magneton, S is the 
total spin angular momentum operator of the nucleus. In 
this notation, the nuclcon magnetic-moment operator is 
written as p-{N) — /io[Ms(^)S + fiv{N)T3S], while the 
isoscalar part of the operator for A is /j,oA*s(A)S. The 
isoscalar baryon magnetic-moment parameters used are 
the experimental value 



Ais(iV)=/ip+A'n- 0.880 



(8) 



for the nucleon, and the theoretical value from the non- 
relativistic quark model 



^ls{A)=^l,{N)/5 = 0.9A 



(9) 



for the A. This quark model actually gives the same value 
for both baryon isoscalar magnetic moments. I prefer to 
use the slightly larger value shown for /is (A) obtained 
from the experimental nucleon isovector magnetic mo- 
ment of flv{N) =£p - firmn = 4.71. 

As is known [E2|, the orbital magnetization term is 
already contained in Jc, and does not have to appear 
explicitly. 



For specific components in the initital and final nuclear 
states, only a few terms are allowed in the multipole sum 
shown in Eq. (^) by the triangle rule for angular mo- 
menta. For the dibaryon components included in our 
calculations, the relative BB orbital angular momenta of 
initial and final nuclear states {Li, Lf) are either (0,2) or 
(2,0). Then only the J = 2 multipole term appears for 
the Coulomb target factor. 

The situation for the transverse RME's is slightly more 
complicated. It is controlled by the spatial RME 



(L' = o|| / dn^[Y,in^). 



\L = 2) 



(10) 



•; = 1.3 



where p = p + (q/4) is the mean of the initial relative 
BB momentum p and its final value p + (q/2) after the 
absorption of a virtual photon of momentum q. Part of 
the photon momentum, q/2, goes into the recoil of the 
dibaryon. The functions fn{q) are not needed at this 
point, and will not be given. The important features are 
that the operator must be a spatial quadrupole operator, 
and that k, can only be 1 or 3. 

Eq. ( p^ has the consequence that the only operators 
contributing to the transvers RME for the dibaryon wave 
functions used in the present study are (a) the convective 
and orbital magnetization current terms in (b) the 

spin magnetization current term in Tj\ with J ~2 and 
3, and (c) the spin magnetization current term in T^^^. 



III. CONTRIBUTIONS OF THE DEUTERON 
AA(^Di) COMPONENT 

The d* dibaryon will be treated in this paper as a pure 
AA(^53) state, using the two-centered Gaussian wave 
function parametrized in ||l5|| as a sum of Gaussians. This 
wave function will be specified with other wave functions 
used in this paper in the Appendix. 

The inelastic production of d* from the deuteron 
AA(^_Di) component in the initial state by elastic eA 
scattering is described by diagram |l|a. The Ar gonne- 
28 (A28) deuteron AA(^Di) wave functions |23f| are 
used, expanded in harmonic-oscillator wave functions 
(HOWF). 

However, it is sufficient to give explicit expressions for 
the RME's appearing in the differential production cross 
section only for single-term wave functions such as the 
single Gaussian 



V'd* (p) 



7V*e 



for d* . Here 



N* = {n(3*') 



*2\-3/4 



(11) 



(12) 
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For the deuteron AA(^_Di) component, a one-term form 
of the HOWF is 



Mp, AA/ Di) = /fWiV2e-^ 3^2m(p) , (13) 

where Paa7 is the AA(^_Di)-statc probabihty of the 
deuteron, 



No = 



16n 

15^ 



1/2 



2^-3/4 



(14) 



and 3^2m(p) — p^Y2m{p) is a sohd spherical harmonic. 
Then the Coulomb RME is 

{d*\\M2{q)\\dr) = y^{d*,S\\M2{q)\\dj,D) 

= -^iXqr.fiq), (15) 
where dr stands for the deuteron AA(^£'i) component; 



A = 



2{a*+b) ' 



f{q) - VP^N*N2 



a* +b 



3/2 



(2/3 



(2/3' 



c = A6 + ao; ao = r^/S = 0.12/r 



(16) 



(17) 



(18) 



The ao term takes care of the baryon form factor at the 7- 
baryon vertex. If the A is assumed for simphcity to have 
the same size as the nucleon, the same Gaussian form 
factor appears in all terms of the production amplitude 
in the nonrelativistic quark model. Additional baryon 
form factors that should appear are already included in 
the BB interactions themselves. 

The transverse RMS's can be separated into convcc- 
tive, and spin magnetization terms: 



(^11^02(9)11^7) 



V2 1 



47r Mb " 
-(Ats(S)A*oA)/(<z) 

x[5'A(g) -5.^(9)] , 



(19) 



Jl7| , p3t , denoted below as potential A28. The approxi- 
mate wave function used is given in the Appendix. 

It is convenient to present the calculated results not as 
angular distributions which contain strong dependence 
on the electron energy, but as the effective T factor 

Teff = Tff -f {vTlvL){Ti + Tr^) ■ (22) 

The Coulomb T factor is a function only of the three- 
momentum transfer q, and is energy-independent in our 
model. The transverse terms contain the kinetmatical 
factor vt |vh^ but the energy dependence from the angle- 
dependent term is small except at low electron energies. 

The weak energy dependence of these effective T fac- 
tors is explicitly illustrated in Fig. ^ for a d* mass m* 
of 2100 MeV and for three different electron energies. 
The A28 d(AA''i:»i) wave function is used. The three- 
momentum transfer has the energy-independent minimal 
value of (7niin ~ 238 MeV/c at this value of m* . The mo- 
mentum range shown corresponds to an angular range of 
about 90° (15°, 2°) for 1 GeV (4 GeV, 27.5 GeV) elec- 
trons. 

In the transverse electric T factor, the amplitude from 
the magnetization current is a factor of 2.5 or 3 larger 
than that from the convective current. This means that 
the electric T factor is larger by roughly an order of 
magnitude when the contribution from the magnetiza- 
tion current is included. 

On the other hand, the transverse magnetic T factor is 
very small at small angles, but has become 1/3 as large 
as the transverse electric T factor at 9 ~ 1 GeV/c. 

The calculated cross section for diagram Ha is propor- 
tional to the probability Paa7 of the AA('Di) compo- 
nent of the deuteron, other things being equal. Besides 
the A28 model, the Argonne group has constructed a 
weaker model with Paa7 = 0.23%. The coupled-channel 
models E and F constructed by Dymarz and Khanna p^ ] 
have Paa7 « 0.1 - 0.4%, and a total Paa « 0.4 - 0.5%. 
The relativistic field-theory model studied by Ivanov et 
al. H gives a total Paa of only 0.08%. 

These theoretical estimates are consistent with the best 
experimental information on the total Paa , namely that 
it does not exceed 0.4% at 90% CL from the null result 
of a bubble chamber search for the spectator A++ in a 
vd knockout reaction 



where P = A in Mb and /is(P), and 

5A(g) = ^^|^, gM^>^^^<f. ^ = ^-A; (20) 



(d*||f",r(g)||dr) = (d*||f^r('?)l|rf7)=0. 



(rf*rsT'(9)IM7) 



: Mo 



/8/7g^,(P)5.(g)/(g), 



(21) 



I use a 3-term approximation to the deuteron AA(^Z?i) 
wave function of the Argonne 28-channel potential 



IV. INCLUSION OF THE NN(^Di) COMPONENT 

OF D* 

Although diagram ^ has a structure rather similar to 
that for diagram |l]b, its production amplitudes are much 
harder to calculate because of three complications: (1) 
The wave function in the minor component d* {NN ^ D3) 
needed in the calculation is not to our best knowledge 
readily available in the literature. It has to be evalu- 
ated ab initio. (2) This wave function is very sensitive 
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to the spreading width T of d*{AA^S3) into NN chan- 
nels. Although the production cross section turns out 
to be only mildly dependent on F, it is desirable to use 
the width consistent with the assumed input dynamics 
in the calculation. This too has to be evaluated. (3) The 
wave function has sharp kinks near pq where the nucleon 
energy En{pq) has the value m*/2. We find it simpler 
to perform a one-dimensional integration for the produc- 
tion amplitudes instead of using the harmonic-oscillator 
expansion described in the last section. 

The calculation is thus broken up into three major 
steps in order to handle the three complications listed 
above. 



A. Perturbative treatment of wave-function 
components 

Our first concern is to estimate the accuracy of a per- 
turbative treatment of d*{NN^D3). This is done by 
first examining a similar perturbative generation of the 
d(AA^_Di) from d{NN^Si) using the Argonne-28 po- 
tential from [0: 



,PiAA'D,,p) = {AA'Di,p\d) 
{AA' D,,p\^V\d{NN^S^)) 



1 



{AA''Di\Ois{ci)\NN^Si) 



md-2EA{p) {2TTf 

xvis{q){p + (l\d{NN''Si))d^(l. (23) 

Here AH — rud — Hq, with Hq the unperturbed Hamil- 
tonian and the deuteron mass, Ea{p) is the nonrela- 
tivistic A energy, C'i8(q) is the NN ^ A A tensor- force 
operator 5}"(q)Ti • T2 in the A28 notation with the 
N A isospin operator. The subscript 18 refers to 
channel 18 of A28. The momentum-space potential is 



visiq) 



,.ni 



47r 



^'tO . 3 „2 j_ ,,2 



fi'^ q-^ + /I' 



f{q). 



(24) 



where ^ = 138 MeV is the pion mass, v^^l = 14.91 MeV, 
and f{q) is a numerically generated cutoff function that 
decreases from 1 at g = to at g = 00. 

Equation (|2^) can be simplified to a one-dimensional 
integral if the deuteron NNi^Si) wave function is ex- 
panded as a sum of Gaussians. It is sufficient to give the 
final expression for a single Gaussian wave function of 
the form of Eq. ([ll]), but with falloff parameter /3 and 
an associated normalization constant A^o if the S-state 
probability were 100%. Then 



i,j{p) = ^{AA^Di,p) 



Hp) 



rrid - 2Ea [p) 
where the subscript 7 refers to the spin degeneracy. 



(25) 



h{p) = ^Noe-^'/'^'J v,siq)e-^'^''' 

xj2{ix)q^dq, (26) 

and X = pqj 0^ . 

The actual calculation is made with a three-term fit to 
the A28 deuteron NN{^Si) wave function given 
in the Appendix. The resulting deuteron AA(^Z)i) wave 
function calculated from Eq. (|2^) (long dashed curve) is 
compared in Fig. ^ with the actual A28 wav e fu nction 
(solid curve) in the fitted form given by Eqs. ( |A1| )-(|A^). 
The wave function is so defined that J tlj^{p)p^dp gives 
the fractional normalization P7 in this AA(^i3i) state. 
The calculated value in the perturbative approximation 
is P7 = 0.65% compared to the exact value of 0.419%. 

One can see from Fig. ^ that the inaccuracy in the per- 
turbative treatment comes from the overestimate of the 
wave function at all momenta larger than « 1.3/Lt, where 
fi = 138.0 MeV is the average pion mass used in these 
calculations. The error seems to arise primarily from 
the neglect of a short-range repulsive central potential in 
the AA(^ Di) channel. Such a repulsive potential would 
have reduced the wave function at small inter-baryon dis- 
tances. The size of the error in our perturbative treat- 
ment in this "^Di state seems to be much greater than the 
results found by [|^. We shall see that there are much 
greater uncertainties elsewhere in the present calculation. 
Consequently, I shall consider the accuracy of the per- 
turbative treatment adequate for the present qualitative 
study. 

The same perturbative method is now used to obtain 
the NN{'^D3) component of d*. An added complication 
appears because the AAC^^s) component of d* is a bound 
state embedded in the A^A^ continuum. The calculation 
will require the use of a spreading width into the A^A^ 
continuum (which is essentially just the total width F of 
d*) and a principal- value Green's function. The matrix 
element of the tensor-force operator is also different from 
Eq. (E3). The final result is rather similar to Eq. (Eq): 



Mp) = {NN''D3,p\d* 



AE 



{AEY + (F/2): 



:h*{p). (27) 



where AE{p) — m* — 2Em{p), 



xj2{ix*)q^dq, (28) 

and X* = pq/P*'^. 

In Eq. (|^, a rapidly changing factor AE/{{AE)'^ + 
(F/2)^) has been separated from a function h*{p) that 
does not depend explicitly on m* or F. AE changes sign 
as the nucleon momentum p increases above that value 
po at which AE{p) vanishes. Above po, AE is nega- 
tive so that the perturbative wave function "ipaip) tends 
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to have the same phase relation to the driving AA(^5'3) 
wave function as tpr^p) is to its own driving wave func- 
tion. This means that diagrams |l|b and c tend to interfere 
constructively at large momentum transfers, but destruc- 
tively at small momentum transfers. 

The falloff parameter f3* and normalization constant 
Nq appearing in Eq. are those for the single- 

Gaussian approximation to the d*{AA'^S3) wave func- 
tion. 

In the actual calculation, a three-term Gaussian fit to 
the two-centered d* wave function is used. The results 
for potential A28 are shown as long dashed curves in Fig. 
^ for h* (p) and the radial part V'3 (p) of Eq. (^) calcu- 
lated with m* = 2100 MeV and T = 100 MeV. This 
spin triplet wave-function component has the perturba- 
tive normalization of P3 — 0.038, or 3.8%. This normal- 
ization is very sensitive to the choice of the decay width, 
increasing to 0.55 if the width is decreased to 10 MeV. 
The rapid increase comes from the roughly 1 /F behavior 
of the amplitude of the kink where the wave function os- 
cillates rapidly from a positive value to a negative value 
as the momentum p goes through the zero of AE. 

The great sensitivity of the "^D^ wave function to the 
width F does not mean that the inelastic production am- 
plitude increases just as dramatically. The wave function 
changes sign in the kink, leading to much cancellation 
in its contribution to the inelastic RME. Consequently, 
most RME's increase by only 25% when the width de- 
creases from 100 MeV to 10 MeV, and approach stable 
values for smaller widths. 

Although the dependence on F as described above is 
technically correct in a narrow sense, it is also counter- 
intuitive in that a smaller width with weaker coupling 
to the NN channel somehow leads to a stronger ^D^ 
component. The truth is that in the discussion just given, 
the width is treated as if it were a free parameter when 
it actually is not. Rather once the dynamics is chosen, a 
certain width is implied, and must be used in Eqs. (|25| ) 
and dH). 

Hence it is desirable to have a more consistent treat- 
ment. I shall first calculate the width for each input 
potential at each dibaryon mass m*, and then use this 
calculated width to calculate admixed components. 

B. Decay width of d* 

The decay width of d* into the NN {ttNN, dnn) chan- 
nel has been estimated in ( ||l^, |^). It is clear that 
the NN channel is by far its dominant decay channel; it 
is the only channel that has to be included in the present 
study. The d* decay widths into NN for the A28 po- 
tential ijl^ and five different Bonn potentials |^,|l8| are 
easily calculated by the same method. The results are 
shown in Table 1 for m* — 2100 MeV and a single Gaus- 
sian wave function for d* of radius r* = 0.7 fm. 



The potentials are ordered in Table 1 in decreasing 
values of their deuteron D-state probability Pd- The 
differences in Pd have been obtained by adjusting the 
amount of short-range repulsion in the tensor potential 
by changes in the coupling constants and baryon form 
factors. 

These potentials differ from one another in other in- 
teresting ways. Potential A28 has a NN AA ten- 
sor potential from tt exchange with an overall strength 
f^AN/ f^NN relative to the NN tensor potential given 
the Chew-Low value of 4 |2^ . The potential function has 
the Yukawa form corresponding to a meson propagator 
l/uj'^{q) in momentum space, where uj^{q) = 
Its short-distance cutoff is specified in coordinate space 
using a function that bends the potential back to the 
origin to simulate the contribution of p exchange ||2^ . 

The Bonn potentials use nuclcon form factors for the 
cutoff, with a monopole form for the tt vertex and usually 
a dipole form for the p vertex. The cutoff masses 
used are shown in the table. The meson (m) coupling 
constants used here for the Bonn-A, B, and C potentials 
jlSf , denoted here as potentials BA, BB and BC, and for 
the "Full" Bonn potential FB are the quark- model 
values based on the ratio fmAN/fmNN = 72/25 [ p8| . 
Potentials BA, BB and BC are relativistic momentum- 
space potnetials. FB is a relativistic model with energy- 
dependent meson propagators. The coupled-channel III, 
denoted here as CC3, is nonrelativistic and uses a strong- 
coupling ratio of about 4.4. 

The second numerical entry given in the first line 
for each potential is the width calculated with the tt- 
exchange only of the potential. The result is clearly cor- 
related with cutoff mass A^ in the Bonn potentials. It 
decreases roughly monotonically from 130 MeV to only 
30 MeV as we go from potential BC to BB, BA, Full 
Bonn (FB) and CCS potentials as A^ decreases from 3.0 
to 0.8 GeV. 

When the p-exchange contribution is added, the re- 
sulting width, shown as the third numerical entry in the 
first line, decreases dramatically to the range of 57 to 
20 MeV. The trend appears to be correlated with the 
decreasing value of the deuteron D-state probability Pd , 
except for the special case of potential FB discussed sep- 
arately. It is of course the differences in Pd caused by 
different partial cancallations by the p-exchange poten- 
tial that distinguish between these potentials in the first 
place. 

The result for the FB potential is different because 
the NN AA potential generated by a simple quark- 
model prescription is not that shown in the first line of 
the results given in the table for the FB potential, but on 
the fourth line, where the width of 17 MeV is in rough 
agreement with the value for potential BA and with the 
trend dictated by Pd- Note that the pAN form factor 
is a monopole, a difference noted in the table by the 
entry [(n^,np) —] (1,1) under the column "Specials". 
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The results of line 4 for the FB potential are also the 
results reported previously in 

In the actual FB potential, the meson-AA^ vertices are 
further modified: First, the pAN form factor is changed 
from a monopole to a dipole form with no change in the 
cutoff mass Ap. This change of the form factor causes the 
p contribtuion to be reduced significantly, thus causing 
the "unexpected" increase in the decay width. 

A second change made in the actual FB potential is to 
decrease the cutoff mass A^ from 1.3 to 1.2 GeV. This 
change is dictated by the need to prevent the uncorre- 
lated 27r contributions to the potential in the lower partial 
waves from becoming unmanageably large ||2^ . a situa- 
tion that is not yet encountered in our lowest-order cal- 
culation. As far as our lowest-order result is concerned. 
Table 1 shows that the effect of this change is quite mi- 
nor. 

The widths for potential CCS also stand out from the 
trend. Comparing them with those for potential BB, 
which has a comparable Pd, we can see that particularly 
the TT contribution has been reduced significantly by the 
use of the much smaller cutoff mass, thus giving an abnor- 
mally small width in perturbation theory. Presumably, 
the channel- coupling potentials in CCS will give rise to 
larger contributions from the long-range 7r-exchange po- 
tentials in higher orders. Hence it is likely that the large 
difference seen in the lowest-order calculation reported 
here will not persist so noticeably in higher orders. 

To summarize, the d* decay width appears to be con- 
trolled to a good extent by Pd, and has a value in the 
range 15-60 MeV. In comparison, the value of 100 MeV 
for potential A28 is more appropriate to 7r-exchange only. 
In other words, the simulation of p exchange contained in 
it has not been very effectively in the long-distance part 
of the tensor potential that controls the decay width. 

I now turn to the many remaining entries in Table 
1. The second line for each potential, or for each group 
within a potential, contains results obtained by replac- 
ing each cutoff factor 1/[1 + (g^/A^)] by the Gaussain 
exp[— (g^/A^)]. This change, marked in the column "Spe- 
cials" by the symbol G, improves the high-momentum be- 
havior of various momentum integrals and also reduces 
the overestimate of the high-momentum wave function 
in our perturbative treatment. Its influence on our decay 
widths is minimal, except in potential CCS where the 
small cutoff mass A^r allows a relatively large effect. 

In using the potential models discussed here, we either 
use the given coordinate-space potential v{r) or the stan- 
dard meson propagator l/a;^((7), plus corrections from 
vertex form factors. The exception is the coupled-channel 
model CCS of ||l^, for which we use the off-shell propa- 
gator 

V{q) ^ u;-^q)[u;{q) + AM]-\ (29) 
defined in the potential. Here AAI = AIa — Mm- 



The question could be raised as to whether one should 
use the off-shell propagator in these cross-channel poten- 
tial as a general policy. Taken by itself, the answer is 
probably no. The reason is that the standard propaga- 
tor leads to energy-independent Yukawa potentials that 
are actually superior to energy-dependent potentials aris- 
ing from the use of off-shell propagators. It is not obvi- 
ous, but the use of an energy-independent potential ac- 
tually corresponds to the inclusion of that component of 
the two-baryon state having a virtual meson "in the air" 

Although the use of the off-shell propagators is not rec- 
ommended, it is clear that their use will reduce the de- 
cay width calculated with the standard propagators. To 
study the size of the reduction, I now repeat the calcu- 
lation with off-shell propagators. The results are shown 
on the third line of each group, and marked by the sym- 
bol AM under the "Specials" column. We see that this 
off-shell effect is quite large, especially on the 7r-exchange 
contributions. As a result, there is much closer cancella- 
tion than before between the tt and p contributions, thus 
giving rise to a much smaller net decay width. 

In subsequent calculations in this paper, I shall use the 
A28 potentials to give an extreme case, and the BB/G 
potential for a more representative example, of the results 
expected for the electroproduction T factors. The decay 
widths calculated for these two potentials are shown in 
Fig. H as functions of the dibaryon mass m* . Gaussian 
vertex form factors are used in potential BB/G. 

With these calculated widths, the '^£'3 component can 
now be estimated more reliably by perturbation theory. 
The results are given in Fig. ^ as long dashed (solid) 
curves for potential A28 (BB/G) using m* = 2100 MeV 
and the calculated value of F = 100 MeV (S8 MeV). Both 
the smooth part h* (p) and the complete radial wave func- 
tion ^{NN ^D:},p) of Eq. ( p7|) are shown. The pertur- 
bative normalization in this state for potential BB/G is 
P3 = 4.5%. 

The perturbative result for 'i(j{AA'^ Di,p) of Eq. (??) 
calculated for potential BB/G with to* = 2100 MeV and 
F = S8 MeV are also given in Fig. ^ as a thick dashed 
curve when only the 7r-exchange part of the potential is 
included in the calculation, and as a thick solid curve 
when the p-exchange contribution is also included. The 
large reduction in the wave function caused by the inclu- 
sion of p-exchange is worthy of note. The perturbative 
normalization in this state for the complete BB/G po- 
tential is P7 = G.S8%. 

C. Electroproduction T factors 

The inelastic production amplitudes from the dom- 
inant deuteron NN{^Si) component to the NN{'^D3) 
component of d* , as described by diagram ^ can finally 
be calculated. The major difference from the procedure 
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described in Sec. Ill is that the radial wave function 
ip3 = d* {NN ^Ds) is now too complicated to be expanded 
readily in terms of harmonic-oscillator wave functions. It 
is kept in numerical form so that the final result, instead 
of being entirely analytic, now requires a one-dimensional 
numerical integration. 

To simplify the calculation, the deuteron d{'^Si) wave 
function is expressed as a sum of three Gaussians. How- 
ever, results need to be given only for a single Gaussian. 
For this purpose, I use the notation of Eq. (pi]), but with 
un-starred parameters. All production amplitudes then 
have the generic form 

{dCSi)\\d{q)\\dl) - co7VoC-('/')'/''''e-"««'/2 

X / e-P'/''''Mp)foip,q)p'dp, (30) 



but different integrand functions 



fo{p,q)^j2, forMz; 
_ V6 2/3^ . 
" Mn q 
^ ^is{N) ci 

^ fi,{N) 2 
Mn 



4^ 

q 

10/32 



J2 - IPJI 



32 - IPJI 



(31) 



Here ji = ji{ipq/ 0^) is a spherical Bessel function, and 
Co = v^, Ci = 1. 

The production amplitudes from diagram la have a 
similar form, but with cq = VS, ci = 4, and of course 
different radial wave functions. The baryon magnetic- 
moment parameter that appears is now /^s(A) instead 
of fis{N). I have verified explicitly that they give the 
same numerical values as the oscillator expressions given 
in Sec. HI. The difference in the numerical coefficients Cq, 
comes from the recoupling of angular momentum. These 
alone favor the amplitudes form diagram ^3 by a factor 
of roughly \flji^ which is a ratio of two 6j-symbols. 

To perform the calculation for potential BB/G, I use a 
three-term approximation to the BB deuteron NNi^Si) 
wave function [p^ . The resulting perturbative deuteron 
AA(^Z?i) wave function calculated from Eq. (25) is fit- 
ted to a 4-term harmonic-oscillator form. Both fitted 
wave functions are given in the Appendix. This fitted 
AA(^Z?i) wave function is the one shown as a thick solid 
curve in Fig. |[ 

The calculated effective T factors are shown in Fig. ^ 
with the curves defined in the legion. 

we can see that electroproduction RMS's from diagram 
|l|a (in the case of A28) or b (BB/G) and diagram |l]c inter- 
fere destructively at small q values. At g = 300 MeV/c, 
for example, the production amplitude that contributes 



to the Coulomb factor from diagram |l]; for potential 
A28 is opposite in sign and about 3.6 times larger than 
that from diagram |l]a. Of the increase, a factor \fTJ% 
comes from a 6j-symbol, while the remaining factor of 
2.4 comes from the wave functions in a radial integral. 
Hence is about 7 times that for diagram |l|a alone. 

The production amplitude from diagram |l]c varies 
more strongly with the momentum transfer because 
of the much stronger momentum dependence of the 
d*{NN^D3) wave function it contains. This produc- 
tion amplitude changes sign just below q = 1000 MeV/c 
where the curve for diagrams [^a and c crosses that for 
diagram |l|a alone the second time. This is above the in- 
terference zero that can be seen in Fig. |. After the 
sign change, the amplitudes from the two diagrams add 
constructively. 

In the effective transverse T- factors for potential A28, 
the production amplitudes from the two diagrams are of 
comparable magnitudes but opposite signs at low mo- 
mentum transfers. The resulting destructive interference 
is very severe, causing the total contribution to Tcff to 
be some two orders of magnitude smaller than that from 
diagram [^a alone. The interference eventually becomes 
constructive at roughly the same momentum transfer as 
for the Coulomb factor. 

The results for the BB/G potential is qualitatively sim- 
ilar, but the Coulomb factor is smaller at the smaller 
momentum transfers. 



V. RESULTS AND DISCUSSIONS 

For all electron energies greater than a few GeV, the ef- 
fective T factors are essentially energy independent. The 
remaining factors in the integrated differential produc- 
tion cross section da/dQ of Eq. (|l]) do depend on the 
energy, however. The differential cross section then in- 
creases roughly as the square of the scattered electron en- 
ergy at comparable four-momentum transfers and small 
lab angles. This means that the cross section tends to 
increase quadratically with the incident electron energy 
as the latter increases. 

The resulting angular distributions are shown in Fig. 
1^ for 1 GeV electrons. The result for potential A28 is for 
diagrams |l| a and c, while those for potential BB/G are 
for diagrams ^b and c. The interference between the two 
diagrams can be visualized more readily by also exam- 
ining Fig. ^ which show additional results for potential 
BB/G from the diagram |l]c alone. 

For the dibaryon mass m* — 2.1 GeV, the sharp local 
minima near 60° comes from the amplitude zero in the 
Coulomb factor. When diagram ^ appears alone, 
this zero comes from the sign change in the d*{NN^D3) 
wave function. When both diagrams 0b and s are in- 
cluded, it comes from an interference zero in the sum 
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of their amplitudes. The two diagrams interefere de- 
structively below the sharp minimum, and constructively 
above it. For example, the integrated cross section at the 
lab angle of 30° is 0.05 nb/sr for potential BB/G, but 0.11 
nb/sr for diagram |l|c alone. 

The over-estimate of the cross section by potential A28 
can be seen in Fig. |^. At the lab angle of 30°, the cross 
section of 0.09 nb/sr for potential A28 is about twice as 
large as the value for potential BB/G. Other than this, 
both cross sections show the same decreasing trend with 
increasing lab angle. 

The Mott cross section ctm and the target recoil factor 
/rec do not depend on the dibaryon mass m*. However, 
the magnitude of four-momentum transfer decreases as 
m* increases, while that of the three-momentum transfer 
increases. This leads to a rapid decrease in the kineti- 
matical factor vi^. 

The dynamical factor Toff also changes with m* . This 
is partly because of the three-momentum transfer ap- 
pearing in it and partly because the d* wave function 
d* {NN ^ D3) itself also changes. The interference be- 
tween the two diagrams becomes rather complex. The 
Coulomb factor decreases rapidly with increasing m* . 
At m* — 2.3 GeV, the minimum in the angular distribu- 
tion at 135° is caused by the interference zero in both 
and T^. At to* = 2.4 GeV, the shallow minimum 
near 90°, comes from the combined effects of an interfer- 
ence zero in T'^ at 100° and an interference zero in at 
50°. The destructive interference between the two pro- 
duction amplitudes at small momentum transfers tends 
to become more severe as to* increases. 

Fig. |7| shows how rapidly the cross section decreases 
with increasing to*. At the lab angle of 30°, the inte- 
grated cross section with both diagrams present, which 
at 0.05 nb/sr is already small for m* = 2.1 GeV, now 
falls down to only 0.0005 (0.00006) nb/sr for to* = 2.3 
(2.4) GeV, as shown in Fig. 0. 

The very small calculated cross sections obtained for 
both larger angles and larger to* suggest that the present 
lowest-order picture might not be adequate under these 
circumstances. A much more elaborate calculation that 
includes higher-order diagrams as well as additional 
production processes important at these larger three- 
momentum transfers might have to be considered. 

The result reported here can be compared with a recent 
calculation by QSW jl^ , who include diagram ^ and a 
production amplitude going to the iVA^*(1520) compo- 
nent of d* . They find that diagram |l|c dominates the 
calculated cross section for 1 GeV electrons scattered to 
a lab angle of 30°, and obtain a preliminary value for the 
integrated cross section greater than 10 nb/sr there. 

In order to compare with QSW, we show in Fig. ^ our 
result for diagram Ic alone at the same m* — 2.1 GeV 
but calculated with potential BB/G. Our cross section is 
only 0.11 nb/sr at 30°, two orders of magnitude smaller. 

One difference between the two calculations is that the 



baryon-baryon interaction used by QSW comes from |Tc| ] 
and contains only the 7r-exchange contribution. In poten- 
tial A28 and certainly in the Bonn potentials used here, 
the baryon-baryon tensor potentials are significantly re- 
duced by cancellation against additional p-exchange con- 
tributions. 

To isolate the contribution of 7r-exchange alone, we re- 
peat the calculation for diagram Ic alone with only one 
modification — using the interaction of jlOj (called here 
potential GO) to generate the perturbative d* {NN ^ D3) 
wave function. The potential GO has the coupling con- 
stant = 0.36 (= /*V47r in the notation of |0|), 
roughly consistent with those shown in Table |[ A 
monopole form factor for the ttAA'^ vertex is used with the 
cutoff parameter A^r = 1250 MeV. The ^Si wave func- 
tion used in this calculation remains that of the Bonn 
B deuteron. This is not very different from the Bonn C 
wave function used by QSW. 

The resulting cross section, for to* — 2.1 GeV, is shown 
in Fig. ^ as a dotted curve. The calculated cross section 
at 30° is 0.41 nb/sr, almost four times bigger than the 
BB/G result. This shows that the neglected p-exchange 
contributions do reduce the cross section significantly. 

The situation remains qualitatively the same when di- 
agram la is also included: The cross section for potential 
GO at 30° is 0.29 nb/sr, a factor of 6 greater than the 
value of 0.05 nb/sr for potential BB/G. 

The cross section for diagram Ic alone for potential 
BB/G calculated with m* — 2.3 and 2.4 GeV are also 
given in the figure and compared with the results for 
TO* = 2.1 GeV with only diagram Ic and with both dia- 
grams present. 

After accounting for the difference in dynamical inputs, 
I find a remaining discrepancy of more than an order 
of magnitude between my calculation and that of QSW. 
This remaining difference must be due to differences in 
the calculational methods used. The most important dif- 
ference seems to be my use of a principal-value Green's 
function in my d*{NN^D3) wave function compared to 
the use of an outgoing- wave boundary condition in QSW. 
The QSW calculation thus gives an additional term in the 
production amplitude from diagram ^ that contains an 
energy-conserving (S-function. This additional term has a 
part that is the decay amplitude of c?*(AA into the NN 
channel. I do not include this part as the physical elec- 
troproduction production amplitude. Additional studies 
must be made to understand if this could account for all 
the remaining disagreement. 

Returning to the general problem of calculating the 
electroproduction cross section of d*, there are of course 
the additional uncertainties in the predicted mass and 
structure of d* itself. Given the experimental non- 
observation of 1=0 dibaryons in the mass range 2.00- 
2.23 GeV |l^, future experimental searches and theoret- 
ical studies should probably move to higher masses. For 
electroproduction, the rapid decrease of the cross section 
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with increasing d* mass is a cause of concern if it persists 
in higher-order calculations. 

The theoretical uncertainty in the d* structure remains 
a major obstacle. The didelta model of d* used here is 
very crude, and may not be adequate at short distances 
where the baryons overlap. Other exotic components in 
the deuteron wave function must also be taken into con- 
sideration, especially at the larger momentum transfers. 
The possibility of quark delocalization and color screen- 
ing remains an open question, but the failure to see any 
d* dibaryon in the experiment of is discouraging. 

All these issues make it clear that much more work re- 
mains to be done before a quantitative description of the 
d* production cross section can be achieved, given any 
theoretical model for d* . The preliminary study given in 
this paper does suggest that the electroproduction cross 
sections to d* are likely to be very small. 

Perhaps more interestingly, the existence of d* at low 
masses can now be viewed from a rather elementary per- 
spective. Consider a very naive model of d* contain- 
ing a probability Paa(c^*) of the AA(''S'3) configuration 
and the remaining probability, 1 — Pa Aid*), of a closed- 
channel configuration that does not decay at all into any 
channel. The two pieces of information: 

(a) the best available experimental upper bound of 0.08 
MeV obtained by ||l6| for its decay width into the 
TV TV channel near 2.1 GeV, and 

(b) the best theoretical estimate of its width of about 
40 MeV given in this paper if d* is a didelta, 

if taken literally, would require that PAA{d*) cannot ex- 
ceed something of the order of 0.2%. This seems to 
suggest that any model of d* containing a much larger 
PAA{d*) might already have been excluded by the mea- 
surements of p6| . 

I wish to thank Fan Wang, Terry Goldman, Stan Yen, 
Earle Lomon, Bob Wiringa, Faqir Khanna and R. Dy- 
marz for many stimulating conversations and correspon- 
dence. I am indebted to Stan Yen for asking about the 
resonance contribution to the total cross section measure- 
ment of Lisowski et al. p^ ] that led to the experimental 
upper bound given in the Introduction. 



APPENDIX A: WAVE FUNCTIONS USED IN 
THE CALCULATION 

Many of the wave functions used have been fitted to 
the form 



(Al) 



f3^ = i(3lf]l...). 
The dimensionless expansion coefficients are 
c (C1,C2, ...) . 



(A2) 



(A3) 



In order to emphasize the stronger high-momentum 
components in these short-distance wave functions, the 
range parameters are obtained by minimizing the per- 
centage mean-square deviation. 

The fitted parameters are given in Table ||. Each fitted 
wave function has been renormalized in order to change 
the raw percent probability to the value the original wave 
function has, as indicated in the third line of the table 
for each state. 



where i^iip) is a normalized harmonic-oscillator wave 
functions (HOWF) hke that shown in Eq. (|l3|), with 
the falloff parameter 
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FIG. 1. Electroproduction of d* from the deuteron: (a) 
from the deuteron AA ^_Di state, (b) from the perturbative 
deuteron AA ^ Di wave function, and (c) from the deuteron 
^S'l state to the perturbative d*{NN^D-s) wave function. 



FIG. 2. The effective T factors appearing in the differential 
electroproduction cross section of d* from the deuteron at 
different electron lab energies as calculated from diagram la 
for potential A28. 



FIG. 3. Comparison of the perturbative wave functions 
ip7{p) in momentum space of the deuteron AA'^Di state for 
the potentials A28 and BB/G to the exact value for the A28 
potential. Here fj, = 138.0 MeV is the pion mass. 



FIG. 8. Integrated differential cross sections da/dfl in the 
lab for the production of d* from deuteron for 1 GeV electrons 
calculated for potential BB/G at different dibaryon masses 
m* when only diagram Ic is included. The dotted curve gives 
the result when the perturbative d*{NN^D3) wave function 
is calculated using the potential GO instead. For comparison, 
the result for m* = 2.1 GeV when both diagrams lb and c 
are included is shown as a thick solid curve. 



TABLE I. Decay width r{d*) in MeV in the didelta model 
of the d* for different baryon-baryon interactions. 



Potential 
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FIG. 4. The perturbative wave functions of the 
d*{NN^D3) state for the potentials A28 and BB/G show- 
ing (a) the smooth function h* (p) , and (b) the wave function 
i^aip) in momentum space. Here /i = 138.0 MeV is the pion 
mass. 

FIG. 5. The decay width F of d* as a function of the d* 
mass m*. 

FIG. 6. Comparison of the effective T factors in the electro- 
production cross section of d* from the deuteron as calculated 
from diagrams la and c for potential A28 and from diagrams 
lb and c for potential BB/G at the electron lab energy of 27.5 
GeV with the results for diagram la alone for potential A28. 

FIG. 7. Integrated differential cross sections da/dil in the 
lab for the production of d* from deuteron for 1 GeV electrons 
calculated for potentials A28 and BB/G at different dibaryon 
masses m*. 



TABLE II. Wave functions used in the calculation. 



Wave function 


Parameters 


d*{NN''S3) 


= 2.0753 (1, 1.25, 1.46) fm^^ 




c = (14.8312, -27.4124, 13.3505) 




100.19% 100.00% 


diNN'-'Si) A28 


= 0.04056 (1, 7.25, 169) 




c = (0.5250, 0.5900, -0.0573) 




93.17% ^ 93.31% 


diNN'-'Si) BB 


= 0.02298 (1, 5.15, 21.9) 




c = (0.3192, 0.4859, 0.3486) 




95.04% 95.01% 


d*(AA 'Di) A28 


= 0.1769 (1, 5.0, 15.0) 




c = (0.00118, 0.01537, 0.05407) 




0.4182% -» 0.4190% 


d^AA 'Di) BB/G 


= 0.9395 (1, 3.18, 8.26, 28.0) 


c 


= (0.00316, 0.02245, 0.02422, 0.02955) 




0.379% 0.383% 
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